STARBURSTS, DARK MATTER, AND 
THE EVOLUTION OF DWARF GALAXIES 
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Abstract 

Optical and H I imaging of gas rich dwarfs, both dwarf irregulars (dl) and blue com- 
pact dwarfs (BCD), reveals important clues on how dwarf galaxies evolve and their star 
formation is regulated. Both types usually show evidence for stellar and gaseous disks. 
However, their total mass is dominated by dark matter. Gas rich dwarfs form with a 
range of disk structural properties. These have been arbitrarily separated them into two 
classes on the basis of central surface brightness. Dwarfs with /j-q(B) < 22 mag arcsec -2 
are usually classified as BCDs, while those faintwards of this limit are usually classified 
as dls. Both classes experience bursts of star formation, but with an absolute intensity 
correlated with the disk surface brightness. Even in BCDs the bursts typically represent 
only a modest <J 1 mag enhancement to the B luminosity of the disk. While starbursts 
are observed to power significant galactic winds, the fractional ISM loss remains modest. 
Dark matter halos play an important role in determining dwarf galaxy morphology by 
setting the equilibrium surface brightness of the disk. 



1 Introduction 



Dwarf galaxies come in a variety of morphologies including: (1) dwarf elliptical (dE) - denned 
by smooth elliptical isophotes, and which invariably have a red color; (2) dwarf irregular (dl) 
- characterized by an irregular structure, blue colors and H II regions scattered haphazardly 
over the optical face of the galaxy; and (3) blue compact dwarf (BCD) jOJ - also an irregular 
morphology, differing from dl in that the H II emission is usually highly concentrated into 
one or two high intensity patches near the center. BCD galaxies frequently are given other 
classifications such as Amorphous |35|, |13| (sometimes prefaced with the adjectives "dwarf" or 
"blue"), or Hll galaxies pi[. However, the properties of dwarf {Mb > —18 mag) galaxies 
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having these different monickers are virtually identical [22, 23|, indicating that they represent 
the same physical phenomenon, hence I will refer to them all as BCDs. 

Despite their morphological differences dls, BCDs, and dEs have similar optical structures 



- their radial profiles are exponential, at least at large radii || f|, |51, |3, |22[]. Naturally this 
leads to speculation: are there evolutionary connections between these different morphologies? 
One commonly accepted evolutionary path was expounded by Davies & Phillips J/J. In it, 
the initial morphology is that of a dl. If its ISM manages to concentrate at the center of the 
galaxy a tremendous occurs starburst occurs resulting in a BCD morphology. Such starbursts 
are known to be capable of powering a powerful galactic wind (e.g. || ^1], p5||). If the wind is 
strong enough all of the ISM is expelled resulting in a dE morphology. If some ISM remains, the 
system fades back into a dl, and undergoes a few more dl ^ BCD transitions before eventually 
expelling all of its ISM to become a dE. 

Here I will address the validity of this scenario by piecing together some work that I have 



been involved in [|2|, |23], |26], |3(J (citation implicit throughout this review), as well as 
that of other researchers, that answers some smaller questions. What kind of dwarf galaxies 
host starbursts? What is the effect of starbursts on the ISM of dwarf galaxies? In addition I 
will consider another issue - dark matter (DM). This is the dominant form of mass in dwarf 
galaxies. Are the differences in dwarf galaxy morphology related to DM content? How does 
DM effect the evolution of dwarf galaxies? 

This review is heavily weighted towards gas rich dwarfs and in particular BCDs. In the 
Davies and Phillips scenario they are the active link between dl and dE galaxies, and hence 
should provide the best clues to deciphering dwarf galaxy evolution. The paper is laid out 
as follows: §|2| compares the optical structure of dls and BCDs; §|3] details the H I structure 
and dynamics of two BCDs: NGC 1705 (D = 6.2 Mpc) and NGC 2915 (D = 3.1 Mpc), and 
compares them to dl galaxies; and §^ synthesizes the optical and radio results to form a new 
scenario where DM plays a dominant role in determining the morphology of gas rich dwarfs. 
Throughout this review I adopt distances based on H = 75 kms" 1 Mpc -1 . 



2 The optical structure and classification of gas rich dwarfs 

Figure [l] illustrates typical radial surface brightness (/i) and color profiles of BCDs using 
NGC 1705 and NGC 2915 as examples. The exponential nature of the outer profiles is clearly 
demonstrated by the linear fits in the /i versus radius R plane. While in dwarf galaxies it 
is not clear that the stars responsible for this light are rotationally supported (in the lowest 
luminosity dEs they are not p4j), I will call this structure the disk for brevity's sake and be- 
cause both dls and BCDs contain rotationally supported Hi (§0). The central region usually 
displays relatively blue emission in excess of the extrapolated disk. I will call this structure 
the starburst because it is responsible for the starburst characteristics of BCDs, as is clear 
from the following evidence. Firstly, Ha imaging shows that the most intense H II emission is 
confined to this region in BCDs. Secondly, HST || [T/J and ground based fl], |5j, |26| imaging 



reveal numerous young blue star clusters in these structures. Finally, their blue broad band 
colors indicate they must be due to young stellar populations: they have ages of ~ 10 Myr if 
instantaneous burst models are adopted or ~ 100 Myr if constant star formation rate models 
are adopted. The strong Ha fluxes from the cores are more consistent with the constant star 
formation rate models (because ~ 10 Myr old instantaneous bursts are no longer ionizing). In 
either case this is much less than the Hubble time, thus confirming their starburst nature. In 
comparison, the colors of the disk are typically like those of stellar populations forming con- 



tinuously over a Hubble time (i.e. like dl galaxies, cf. [O, [34]), or a bit redder suggesting an 
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Figure 1: Surface brightness profiles in B and R, with corresponding (B — R) color profiles 
of NGC 1705 and NGC 2915 |§. The dashed lines show exponential fits to the outer 



portions of the /i profiles. Note: the optical profiles of NGC 1705 exclude the light of the bright 
cluster NGC1705-1. 



inactive population. Hence, the disk is best attributed to the progenitor, or host, galaxy. 

Surface brightness profile fitting provides a means to determine both the relative strength 
of the starburst, and the structural properties of the disk. The method is simple. The outer 
portions of the profile are fitted with an exponential. This yields the extrapolated central 
surface brightness /io (which is corrected for inclination), and scale length a -1 of the disk. The 
burst and disk are separated by assuming that the disk remains exponential all the way into 
the center. 

Figure ^| shows the ratio of B band fluxes of BCD bursts relative to their host disks. The 
strongest starbursts are about twice as bright as their hosts. Hence, while starbursts can 
outshine the host disk they are nevertheless modest < 1 mag enhancements to the total B flux 
of BCDs (cf. pm ). Typical flux enhancements are only a few tens of percent. This figure does 
not include upper limits to the burst/disk ratio: about 20% of BCDs have exponential profiles 
all the way into their cores (see below). The mass contribution of the starbursts is even smaller, 
typically < 5%. These are not the ~ 6 mag starburst enhancements proposed to explain the 
excess of faint blue galaxies at moderate redshifts [§]. 

Figure [I] compares the disk parameters /io, of both BCDs and dl galaxies. I reemphasize 
that does not include the contribution of the starburst core. While there is some overlap, 
we see that BCD distribution is offset from that of dls in terms of both /i and or 1 . The offset 
in a~ l is in the sense that BCDs tend to be smaller than dls. However this difference may 
be due to the (loose) selection criterion M B > — 18 mag typically applied to dwarf galaxy 
samples: large scale length, high surface brightness galaxies are not dwarfs. More striking is 
the difference in /z (cf. |32|). Typically /io is 2.5 mag arcsec -2 more intense in BCDs than in 
dls. Structurally BCD disks are very different from those of dls. 

The absence of BCDs on the left half of Fig. |3] is puzzling. Does this mean that dl galaxies 
do not experience starbursts? Examination of the /i and color profiles of Patterson & Thuan 
[ 3~Jfl reveal several dls (e.g. UGC 5706, UGC 7636) with /a profiles of like those in Fig. [I], that 
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Figure 2: (Left) Ratio of B band fluxes of starbursts to host disks in two samples of BCDs 
(symbol correspondence in right panel). 

Figure 3: (Right) Structural parameters of the exponential disks for a variety of samples of 
BCDs, and one dl sample. 



is containing a blue central "high" surface brightness excess relative to the exponential disk. 
This is structural evidence for starbursts in dls. The episodic star-forming nature of dls is 
best demonstrated using color-magnitude diagrams of the nearest ones (e.g. |T4]|). However the 
observed central surface brightness of the bursting dls including the light of this central excess 
is typically fi(B) > 22 magarcsec -1 , much fainter than the central regions of BCDs. While 
dls do experience short duration bursts of star formation, they are pathetic and not usually 
recognized as starbursts because they are not intense enough. 

On a similar note, about 20% of BCDs (data from [IIJ) have \i profiles that are nearly 

exponential and have fairly flat color profiles (examples include Haro 14 [^] and UM483 |4^| ) 
- that is with no structural evidence for a starburst. The colors of these galaxies are consistent 
with fairly long duration ( > 1 Gyr) star formation. These results highlight the importance of 
absolute surface brightness in morphological classification. For example, in the Virgo Cluster 
Atlas [36 galaxies are recognized as dwarfs by their low surface brightness. Similarly BCDs 
are recognized by their high surface brightness. This is largely an unstated, perhaps even 
unconscious, selection criteria. There are clearly low surface brightness dwarfs (e.g. GR8) that 



meet the usual luminosity, emission line, and size criteria for BCD classification ||5[] but are 
invariably classified as dls. From Fig. |3] the BCD/dl dividing line is (io{B) w 22 magarcsec -2 . 
Clearly the presence of a starburst makes it more likely that a dwarf will be classified as a BCD. 
However, such a classification does not guarantee the presence of a strong burst. Nor does the 
dl classification exclude galaxies experiencing strong bursts relative to the host disk brightness. 



3 H I structure and dynamics of BCDs 

Radio array observations of the H I structure and dynamics of BCDs can tell us much about 
star formation feedback (e.g. galactic winds) and allow the distribution of mass (including DM) 
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to be determined. Compared to dls there are not that many H I imaging studies of BCDs; 
in part because it is harder to find BCDs that are resolvable with radio arrays, due to their 
small numbers and usual compact angular sizes. My collaborators and I therefore decided to 
address these issues with Australia Telescope Compact Array observations of NGC 1705 and 
NGC 2915, two of the nearer BCDs. The resultant composite H I and optical images are shown 
in the color section of this volume. Although both galaxies have some kinematic irregularities 
their dominant structures are extended rotating disks which are strongly centrally peaked. 
These are typical properties of BCDs imaged in HI EH, fill E7|. The disk of NGC 2915 is so 



extended that it has the H I appearance of a late type barred spiral. Similar galaxies include 
IC 10 and NGC 4449 



3.1 Starburst - ISM feedback 



Both NGC 1705 and NGC 2915 show evidence of star formation churning up the neutral ISM. 
In NGC 2915, kinks and enhancements in the velocity dispersion map correspond well to Ha 
bubbles and peculiar knots associated with recent star formation. However it does not appear 
that H I is being ejected from the system. In the center of the galaxy, where star formation 
is the most vigorous, am ~ 40 kms -1 which is the same as the one dimensional velocity 
dispersion derived for the DM particles (see Fig. f| below). Hence, star formation appears to 
be maintaining the central H I in virial equilibrium with the DM halo. This suggests that DM 
plays a role in the feedback process: if the starburst energizes H I to have o"hi much larger than 
the halo velocity dispersion, then neutral ISM is thrown into the halo (or beyond) and star 
formation shuts down. 

Jutting out obliquely from NGC 1705's edge-on H I disk out to R ~ 4.5 kpc, is a gaseous spur 
containing 8% of the H I flux (see color plate). Its position angle is similar to that of the optical 
outflow, hence this may be the one-sided blow out of the prominent Ha wind. Its outflow is 
inferred to be predominantly transverse, as is the Ha outflow, hence the expansion velocity V exp , 
or similarly the expansion timescale t exp , is uncertain. The Ha outflow has t exp ~ 10 Myr, like 
the age of the most prominent star cluster NGC1705-1. Adopting this t exp for the Hi outflow 
yields V cxp ~ 300 kms -1 , an average mass loss rate of M. ~ 1.7 MqJt -1 , and total kinetic 
energy E k (gas) ~ 1.5 x 10 55 erg. If t exp ps 100 Myr, closer to the age of the continuosly star 
forming population in the core of NGC 1705, then V exp ~ 30 kms -1 , M. ~ 0.17 M. & yr _1 , and 
-Efc(gas) ~ 1.5 x 10 53 erg. In comparison, the current star formation rate is ~ 0.13 A4q yr _1 . 
Mass loss is at least competitive with star formation in regulating the gas content of NGC 1705. 

Can the energetics of this gas be accounted for by NGC 1705's stellar populations? Using 
the solar metallicity Salpeter IMF (mass range 1 to 100 A4©) population models of Leitherer 
& Heckman |2(J I estimate the mechanical energy output (from stellar winds, supernovae) of 
the young stellar populations in NGC 1705 integrated over t exp = 10 (100) Myr to be E k (*) ~ 
2.3 (11.3) x 10 55 erg. The uncertainties on both Ek(gas) and Ek{*) are probably a factor of few. 
Hence, within the considerable uncertainties, the energetics of the spur are are consistent with 
it being a starburst driven wind. 

The fate of NGC 1705's ejected gas depends critically on V^p. If it is as high as 300 kms -1 , 
then even the large amounts of DM in the system do not gravitationally bind it. If V exp ~ 30 
Myr, some gas will be retained. At least 8% of the neutral ISM has been ejected into the halo of 
NGC 1705 if not out of the system entirely. Hence we are witnessing a significant mass ejection 
event. Nevertheless, even in NGC 1705, a BCD with one of the most spectacular Ha outflows, 
the majority of the ISM is retained in a disk. Even this starburst is incapable of totally blowing 
away the ISM. 
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Figure 4: (Left) Rotation curves of NGC 1705 |p and NGC 2915 [2§ with mass model fits 



The thick solid lines shows the full models, the dotted, dashed, and dot-dashed lines show the 
contributions of the stars, H I disk, and DM halo respectively. The model parameters are listed 
in each panel. The vertical arrows indicate the optical Holmberg radius of each galaxy. 

Figure 5: (Right) DM halo central density po plotted against disk central surface brightness. 
Open squares are from de Blok & McGaugh |§, while diamonds represent NGC 1705 and 
NGC 2915. The top panel shows the results for maximimum disk model fits, while the bottom 
panel shows Bottema disk fits. The circles on the right side of the top panel mark crude 
estimates of po in 12 BCDs with published and unpublished RCs. The dotted line, at bottom, 
is a fit to the Bottema disk results with the relationship log(po) = 0.41og(po) + Constant. 

3.2 Mass distributions 

Figure f| shows the derived rotation curves (RCs) of NGC 1705 and NGC 2915 with mass model 
fits to them. These are the first two BCDs that have been mass modelled. The models consist 
of three components to the mass distribution: (1) the stellar distribution which is given by the 
projected luminosity profile (Fig [I]) scaled by (Ai/ Lb)* - the mass to light ratio of the stars; 
(2) the neutral ISM distribution which is set by the H I profile scaled by a constant 1.33 to 
account for the Helium contribution (i.e. no free parameters); and (3) a dark matter halo. This 
halo is taken to be a pseudo-isothermal sphere with a density distribution given by 

P 1 + {R/R c f 1 ' 

where the free parameters are the central density po and the core radius R c . From these, the 
rotational velocity at large R, Voo, and halo velocity dispersion ctq are given by [|T9"|: 



V£ = AttG Po R 2 c = 4.9< (2) 

The models shown are maximum disk mass models, where (Ai /Lb)* is set by the first points 
in the RC, and then held fixed. Maximum disk models are by nature minimum halo models, 
hence po is a lower limit. We also made minimum disk models ((Ad /Lb)* = 0) and "Bottema 
Disk" models where (At/L B )+ is set by the color of the optical disk H. 
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There are a few things to note about Fig. [|. First, the form of the RCs is like that seen in 
normal disk galaxies, consisting of a rising inner portion and a more or less flat RC thereafter. 
Second, the optical extent of the galaxy, as marked by the arrows is contained within the rising 
portion of the RC. In fact the optical extent marks very well the RC knee (or equivalently R c ). 
The optical portions of dls are also commonly contained within the rising portion of their RCs 
S7| . This is where the rotation is almost solid body, and thus shear is minimized, enhancing the 



ability of clouds to form stars |16] . Thirdly, in both cases DM dominates the mass distribution, 
even within the optical radius of the galaxy. In comparison, the stellar component has a mass 
equal to or less than the neutral gas disk. 

Overall, the global dynamics of BCDs appear to be similar to dls: they are dominated by 
rotating disks with normal looking RCs. A distinction between the two types is seen when the 
DM halo densities po are compared, as shown in Fig. [5]. Central densities found by maximum 
disk and Bottema disk fits are shown in separate panels. The comparison sample is taken 
from de Blok & McGaugh ||, and includes only galaxies with Mb > — 18 mag. It is com- 
prised mostly of dls (left side of panel), but also includes three low luminosity spirals having 
fJ<o(B) « 22magarcsec -2 . This comparison shows that NGC 1705 and NGC 2915 have two of 
the highest po measurements of any dwarf galaxies. In order to check that these galaxies are 
typical, I crudely estimated po from the central velocity gradient for 12 BCDs with published or 
unpublished RCs, and plotted them as circles at arbitrary p in the top panel of Fig. [| These 
estimates are upper limits, since the contribution of the baryonic components to the velocity 
gradients have not been removed. Nevertheless, the comparison indicates that NGC 1705 and 
NGC 2915 have normal po f° r BCDs. Figure [| shows a weak but noticeable correlation be- 
tween log(po) and fio(B), with higher surface brightness disks corresponding to higher po halos. 
This result holds for both maximum disk and Bottema disk solutions, and was first noted (for 
Bottema disks) by de Blok & McGaugh § , who show that it extends to high luminosity disks. 
Here we show that the correlation also includes BCDs. 



4 Evolutionary Connections 

The correlation in Fig. [| can readily be explained by considering the response of a self gravi- 
tating disk immersed in a DM halo core of constant density po, i.e. where the rotation curve is 
linearly rising. Then, the angular frequency Q, p , and dynamical time t^ yn are constant with 
radius and related by 



and Toomre's disk stability parameter Q is given by 



Q - (4) 

where a g is the gas velocity dispersion and E ff is the disk surface density. Lower values of Q 
indicate a higher degree of self gravity. In normal disk galaxies Q typically remains radially 
constant at Q « 2, and star formation correlates with regions of somewhat lower Q |T7| , |T2"| . 
This suggests that star formation regulates disk structure to retain constant Q. Kennicutt [18 
finds that for normal disk galaxies the star formation rate per area is given by 

oc (5) 

^dyn 
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As noted in §|| the optical disks of dwarf galaxies typically have colors indicating a long duration 
of star formation. Let us assume that these are normal disks and hence form stars with the 
above star formation law and at a constant and universal Q. Then equations ||], |] and || can be 
combined to yield 

oc a g p . (6) 

For a constant star formation rate population and a sufficiently blue passband (e.g. B) the 
linearly surface brightness is proportional to E*. Therefore, for DM dominated galaxies we 
expect a simple correlation between linear surface brightness and po- This is consistent with 
the observed correlation, as shown by the dotted line in Fig. ||. Meurer et al. [215 find a similar 



correlation between surface brightness and po holds in the center of starburst galaxies. However 
for them it is normal baryonic matter that dominates po rather than DM. 

In essence, the central mass density determines the equilibrium star formation rate of the 
embedded disk. Following the discussion in §[|, the disk intensity largely determines whether a 
dwarf galaxy is classified as a BCD or dl. As noted earlier, the optical size of dwarfs seems to 
be limited to R c . Hence, both DM halo parameters are important in governing the morphology 
of gas rich dwarfs. Bottema disk mass model decompositions |5] indicate that the two DM 
halo parameters are uncorrelated. This is contrary to cosmological simulations showing DM 
halos forming a one parameter (mass) sequence [31]. The origins and implications of this 
contradiction deserve further investigation. 

Can there be evolution between dl and BCD classes? While some evolution in po may be 
allowed, it is unlikely that there can be enough to change a typical dl into a typical BCD. That 
would require a 2.5 mag arcsec -2 change in p or equivalently a factor of ten change in p . The 
most obvious way to do that is to expand or contract the halo by a factor of e = 10 1 / 3 ~ 2.2 
following a mass loss (wind) or gain (accretion or spin down of an extended disk). The mass 
loss/gain fraction / required to effect this change is given by 

(7) 
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for homologous expansion/contraction that is slow relative to t^ yn [p] ]. A factor of 10 change 
in po thus requires a 55% mass loss or gain. The problem is that there isn't that much baryonic 
mass in a dwarf galaxy. To effect this large of a change would require DM loss or gain. This is 
not feasible if DM is non-dissipative and feels only the force of gravity, as is usually assumed. 
I conclude that there is probably little dl BCD evolution. 

If the ISM were removed from a dl or BCD, it could still plausibly evolve into a dE (or 
dwarf spheroidal) galaxy. However, as noted in § ^ even in a dwarf galaxy undergoing a strong 
starburst with a spectacular galactic wind (NGC 1705), the fractional loss of the ISM is modest. 
If this is typical, it would take on order of 10 bursts to expel all the ISM from a BCD. The 
bursts aren't strong enough, and the ISM distributions are too flattened to allow a single burst 



expulsion of the ISM [pLl ] . Single bursts should have a more profound effect on the chemical 
evolution of dwarf galaxies, since the hot metal enriched ISM is preferentially lost in a starburst 
driven wind fl0 |. This view is consistent with arguments by Skillman & Bender |38| against 



evolution between gas rich and gas poor morphologies occurring commonly at the present 
epoch. The demographics of dwarf galaxy morphologies point to an environmental component 
to their evolution. Gas rich dwarfs are found in low density environments where the frequency 
of external starburst triggers is low. They survive easily. The clock runs faster (more frequent 
triggers) in clusters, and in addition ram pressure stripping would accelerate the removal of gas 
from dwarfs, while tidal truncation of DM halos would assist galactic wind losses. Hence it is 
not surprising that gas poor dEs are found more often in clusters than the field. 
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5 Conclusions 



We are now at a position to re-evaluate the Davies and Phillips [[F] scenario for dwarf galaxy 
evolution. The mechanisms they invoke have clearly been verified. Dwarf galaxies do experience 
starbursts and these can expel some of the ISM. Mass expulsion can rival or surpass lock up 
into stars in regulating the gas content of dwarfs. However, the results of any single burst 
are not so severe. Cataclysmic bursts are not common at the present epoch, and the milder 
bursts that are observed may not be sufficient to change a galaxy's morphological classification. 
The morphology of a dwarf galaxy is largely set by its enveloping dark halo, and is relatively 
impervious to starbursts. 
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